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ABSTRACT 

Many astrophysical problems arising within the context of ultra-high energy cos- 
mic rays, very-high energy gamma rays or neutrinos, require calculation of the 
flux produced by sources tracing the distribution of galaxies in the Universe. We 
discuss a simple weighting scheme, an application of the method introduced by 
Lynden-Bell in 1971, that allows the calculation of the flux sky map directly from 
a flux-limited galaxy catalog without cutting a volume-limited subsample. Using 
this scheme, the galaxy distribution can be modeled up to large scales while rep- 
resenting the distribution in the nearby Universe with maximum accuracy. We 
consider fluctuations in the flux map arising from the finitcness of the galaxy 
sample. We show how these fluctuations are reduced by the weighting scheme 
and discuss how the remaining fluctuations limit the applicability of the method. 

Key words: methods: miscellaneous, catalogues, large-scale structure of Uni- 
verse, galaxies: luminosity function, cosmic rays 



1 INTRODUCTION 

Recent developments in multi- wavelength/multi-messen- 
ger observational techniques often make it desirable to 
calculate the angular distribution of a diffuse flux ex- 
pected from sources with a given spatial distribution. The 
predicted flux distribution may then be used for source 
identification, estimation of the background, etc. Neces- 
sity for such a calculation arises in the context of ultra- 
high energy cosmic rays (UHECRs), neutrino physics, as 
well as gamma-ray astronomy. 

If the sources are extragalactic, their space dis- 
tribution can be derived from the matter distribu- 
tion in the Universe. The latter can be inferr e d from 
galaxy surveys , e.g. lAdelman^ McCarthy et al.l <|2008h : 
ISkrutskie et al.l (|2006l h lJones et al.l (|2005h Good distance 
determination is required to reconstruct the spatial mass 
distribution. Special techniques have been developed to 
minimize the impact of dista nce errors and t o sup press 
the short-scale noise (see, e.g.. lErdogdu et all (|2006l ) and 
references therein). 

The problem of flux calculation has a number of fea- 
tures that make it different from (and easier than) re- 
construction of the full three-dimensional mass distribu- 
tion: (i) only a two-dimensional projection of the three- 
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dimensional distribution is needed; (ii) contributions of 
remote sources are suppressed by the geometrical factor 
r~ 2 and, in many cases, by the flux attenuation due to 
interactions with the ambient matter; (Hi) the smaller 
amplitude of inhomogeneities at larger scales makes the 
contribution of remote sources essentially isotropic; only 
the overall normalization of such an isotropic part has to 
be calculated. These simplifications result in weaker re- 
quirements on the quantity and quality of astronomical 
data in flux calculations, which makes it advantageous to 
by-pass the reconstruction of matter density and calculate 
the flux distribution directly from the galaxy catalogs. Ac- 
curate results may be achieved with substantially smaller 
input. 

Both in the context of mass distribution and in flux 
calculations, a crucial requirement is completeness of the 
underlying galaxy catalog. That is, a volume-limited sam- 
ple is needed which includes all galaxies of a certain kind 
within a given volume. On the contrary, a natural prod- 
uct of an astronomical survey is a flux-limited sample that 
contains all galaxies up to certain apparent magnitude as 
set by the instrumental sensitivity and observation time. 
Volume-limited samples may be obtained from a flux- 
limited sample by cutting away objects that are further 
than a given distance and dimmer than a certain absolute 
magnitude, chosen in such a manner that the resulting 
sample is complete. 

To model adequately the source distribution in the 
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Universe, one requires a galaxy catalog that (i) extends to 
sufficiently large distance (large enough that the Universe 
can be approximated as homogeneous beyond that), and 
(ii) accurately represents the distribution of matter on 
small distance scales. For any single volume-limited sub- 
sample, these requirements work in opposite directions: 
the first calls for a large volume, while the second calls for 
a small volume. To resolve this conflict one could consider 
combining two volume-limited samples, a sparse one to 
cover large distances and a dense one to cover the nearby 
region in more detail. In this paper we discuss a technique 
- termed "sliding-box" technique - which generalizes this 
idea by combining many volume-limited samples. Up to 
a specified limiting distance the whole flux-limited cat- 
alog is used in this process, so that close-by structures 
are mapped out with maximum detail by the dim objects 
in the original catalog. The crux of the construction is 
in an appropriate weighting scheme: We weigh galaxies 
in the flux-limited sample in a distance-dependent way 
so that the progressive incompleteness at large distances 
is compensated by the increasing weight of each galaxy. 
These weights, essentially representing the luminosity per 
galaxy, can be naturally incorporated in flux computation 
algorithms which inherently use some sort of weights to 
account for the fact that remote sources produce less flux 
than nearby ones. 

The sliding-box technique is essentially a method to 
deal with the fact that galaxy catalogs do not contain an 
infinite number of galaxies. The finiteness of a galaxy cat- 
alog unavoidably leads to fluctuations in flux predictions. 
The sliding-box technique strongly reduces these fluctua- 
tions by efficient use of the available data. Nevertheless, 
the remaining fluctuations may still be large enough to 
spoil the accuracy of flux maps modeled from a galaxy 
catalog. We will address this problem in detail and present 
a criterion for the applicability of the sliding-box method. 

An an illustration we will apply the sliding-box tech- 
nique to a su bset of the 2 Micro n All-Sky Redshift Sur- 
vey (2MRS) l|Huchra et all 120091 ). a flux-limited sample 
of galaxies with observed Jf 3 -magnitude m < 11.25 that 
contains measured redshifts for all but a few galaxies Q 
We would like to stress, however, that the technique is 
completely general and can be applied to any flux-limited 
galaxy sample. 

It is worth noting that the issue addressed in this 
paper is related to the well-known astronomical problem 
of reconstructing the luminosity function from a quasar 
or galaxy sampl e that i s limited in ap parent magnitude 
(see, e.g. jFelterJ i| 1977ft ; I Willmerl i| 19971 ) for a comparison 
between different methods and references). The sliding- 
box scheme discussed in this work is an applicatio n of th e 
C~ -method that was pro posed by 
and further developed by IJacksorJ 



Lvnden-Benl (Il97ll) 
1974 ): IChoIoniewskil 



and turtner developed pylJacKsonl l|ltjY4l l; IQhotomewskil 
l|l987t) : lEfstathiou et al.l jl988l ) , to reconstruct the quasar 



1 This subset was kindly provided to us by John Huchra. 
Tailored to model the distribution of galaxies in the field of 
view of a northern-hemisphere cosmic-ray experiment, it does 
not cover the galactic plane with |fe| < 10°, b being the galactic 
latitude, nor the region with 8 < —30°, 8 being the declination 
in J2000 equatorial coordinates. 
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Figure 1. Two nested volume-limited subsamples of a flux- 
limited sample. 



luminosity function. Our formulation of the scheme is tai- 
lored for flux calculations, in keeping with the aim of the 
present study. 

The rest of this paper is organized as follows. In sec- 
tion [2] we discuss the sliding-box technique and present 
an efficient implementation scheme. We also discuss the 
connection between the weights associated with the al- 
gorithm on one hand and the luminosity and selection 
functions on the other. Section [3] is concerned with the 
effect of fluctuations on model fluxes due to the finite size 
of a galaxy sample. In section [4] as an example, we apply 
the sliding-box method to model the flux of UHECR pro- 
tons with energies above 60 EeV from the 2MRS sample. 
We summarize our work in section [5] 



2 THE SLIDING-BOX WEIGHTING 
SCHEME 

2.1 Combining two volume-limited samples 

To illustrate the general idea of the sliding-box technique, 
consider a flux-limited sample that is complete to a given 
apparent magnitude mo. On the M—D plane, where M is 
the absolute magnitude and D is the distance, this sam- 
ple occupies the populated region in Fig. [T] The apparent 
magnitude m of a source is a function of its absolute mag- 
nitude and distance, m = m(M,D). For a given absolute 
magnitude, m increases with distance and reaches the lim- 
iting value mo at a distance D satisfying m(M, D) = mo. 
This determines the line M*(D), the boundary of the pop- 
ulated region in Fig. [1] Beyond this line the objects are 
too dim and the completeness of the sample cannot be 
guaranteed. 

At a given distance Do, only objects with absolute 
magnitude M < Mo ~ M*(Dq) are sufficiently bright 
to be included in the flux-limited catalog. Galaxies that 
are closer than Do and brighter than Mo form a volume- 
limited sample. These are objects in regions marked with 
Si and S2 in Fig. [T] The completeness of this subsample 
follows from the completeness of the original flux-limited 
catalog. 

It is clear from figure [1] that at small distances the 
volume-limited sample Si + S2 contains only a fraction of 
available galaxies which may be insufficient to represent 
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accurately the details of the matter distribution. To im- 
prove the situation, one may construct a denser volume- 
limited sample corresponding to a smaller distance D\ 
(the sample Si + S3 on Fig.[T|. When calculating the flux 
one may use galaxies from S2 + S3 at distances D < Di 
and galaxies from Si at D\ < D < Do. The luminosity of 
a given volume is determined by the number of galaxies in 
the sample S1+S1 contained in that volume. At distances 
D < Di, the same luminosity may be represented in a 
greater detail by galaxies from S2 + S3 provided they are 
assigned smaller "weight" , that is, luminosity per galaxy. 
If the galaxies in the sparse sample have a weight wo 
each, the galaxies in the dense sample S2 + S3 should be 
weighted with 

Si 



^ = s7Ts-3 W0 - (1) 

Here and below we use the same letter to denote the sam- 
ple and the number of galaxies in the sample. At distances 
D < Di, the total weight in the sparse and dense samples 
is the same, SiWq. The difference is that in the dense sam- 
ple it is distributed among a larger number of galaxies, 
and hence the spatial distribution of matter is represented 
more accurately. 

Several volume-limited samples may be combined in 
the same way. In the limit of an infinite number of nested 
volume-limited samples one arrives at the "sliding-box" 
weighting scheme described now. This scheme is essen- 
tially an implemen tation of the C~-method, proposed by 
lLvnden-Bellli|l97ll ). applied to distance D and magnitude 
M. 



2.2 Sliding-box scheme 

Imagine a variable rectangular "sliding box" with one cor- 
ner fixed at zero distance and minimum M (the lower-left 
corner in figure [l]) and the opposite corner moving along 
the line M,(D). At any given position the box defines 
a volume-limited sample. One starts at some maximum 
distance D max ; galaxies that are further than D max are 
disregarded, i.e. assigned a zero weight (this is the part 
of the catalog that is lost). The current weight is set to, 
say, 1. Now the free corner of the sliding box is moved to- 
wards smaller distances. Each time a galaxy exits the box 
through its vertical edge it is assigned the current weight. 
Each time a galaxy enters the box through the horizontal 
edge, the current weight is multiplied by N/ (iV+1), where 
N is the current number of galaxies in the box. When the 
procedure is finished, all the galaxies at D < D max have 
been assigned a weight. 

The main asset of the sliding-box scheme is that the 
weight at a given scale D is computed from a volume- 
limited sample corresponding to distances just slightly 
larger than D; this sample has the maximum available 
number of galaxies and hence the smallest fluctuations. To 
demonstrate the accuracy of the sliding-box method, con- 
sider a direct computational scheme in which the weight 
at distance D is determined from a volume-limited sample 
up to D and a volume-limited sample up to D max (i.e., 
without refining the weights at intermediate distance as is 
done in the sliding-box method). The direct scheme and 
the sliding-box method are equivalent in the limit of infi- 
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Figure 2. Galaxy weights as a function of distance, obtained 
with the sliding-box algorithm and with a direct computation. 



nite galaxies in the original sample. Given a finite num- 
ber of galaxies, however, large statistical fluctuations will 
show up in the direct computational scheme due to the 
sparsity of the volume-limited sample extending to Dmax 
at small distances. This can be seen clearly in figure [2j 
which shows the weights as a function of distance for all 
galaxies in the 2MRS sample computed by the sliding-box 
technique and by direct computation. 



2.3 Relation to luminosity and selection 
functions 

The weight defined by the sliding-box method is re- 
lated to the galaxy luminosity distribution and to the 
selection function characterizing the flux-limited galaxy 
sample. The original construction by iLvnden-Belll (1971) 
was, in fact, aimed at recovering t he quasar luminos- 
ity function (see also IJacksonl 1 19741 ; ICholoniewskil 1 19871 ; 
Efstathi ou et al.lll988l ). To clarify these relations we con- 
sider the problem in general terms. We assume in this 
section that all the samples are very large so that a sta- 
tistical description applies. For simplicity, we also assume 
that the distribution of galaxies in luminosity is space- 
independenlQ, that is the full distribution factorizes into 
a spatial and a luminosity part. The number of galaxies 
with magnitudes between M and M + dM at distances 
between D and D + dD is then expressed as a product of 
two factors, 

dN = \(M)dM ■ v(D)D 2 dD. 

Let N(D, M) denote the total number of galaxies within 
distance D and brighter than M. It factorizes into a prod- 
uct of two cumulative distributions: 



N(D,M) = N(D)L(M) . 



where 



N(D) = f° 
Jo 



dDD u(D) : 



(2a) 
(2b) 



2 The flux non-uniformity arises at z < 1 where one can ne- 
glect the evolution of sources. The contribution from regions 
z ^> I may be non-negligible or even dominant; however, this 
contribution is isotropic. 
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Figure 3. Luminosity distribution \(M) for a volume-limited 
subsample of the 2MRS sample up to 30 Mpc. 



and 



L(M) 



/' 

J — ( 



dMX(M) . 



According to eq. ([T]), the weight function w(D) at 
distance D < D max may be expressed as follows: 



w(D) 



N(D, M*(D ma x)) _ i(M»(D max )) 



N{D,M*(D)) 



L(M,(D)) 



where we have normalized the weights so that w(D n 
1. Equations ([2]) and <[3j) imply that 

/M, (£>) 
\{M')dM' . 
-oo 

This relation can be inverted to read 

w'(D*(M)) dD„ (M) 



\{M) oc 



w(D* (M)) 2 



(2c) 



(3) 



(4) 



(•>) 



where w'(D) — dw/dD and Z)*(M) is the inverse of 
M„(£>) defined by M»(£)»(Af)) = M. Hence we can di- 
rectly infer the galaxy luminosity distribution X(M) from 
the weight function w(D) and vice versa. 

Figure [3] shows the luminosity distribution for a 
volume-limited subsample of the 2MRS sample up to 30 
Mpc derived from the weights, together with the distribu- 
tion directly reconstructed from the data. The agreement 
is excellent. 

Now consider the relation between the weights and 
the selection function (f> sc i(D) in a flux-limited sample. 
The total number of galaxies within distance D in an 
incomplete sample is expressed in terms of the selection 
function as follows, 



n(D) = / dDD 2 v(D)4> sel {D). 



(6) 



Making use of the general expressions in eqs. ((2} one finds 
that 



4>se\(D) = w(D)- 1 • L(M,(Auax)) . 



(7) 



where the last factor is just a normalization constant. The 
selection function derived from the weights as given in 
eq. is shown in figure [4] The model curve is in fair 
agreement with the data, given that it is derived under 
the assumption of a homogeneous Universe. 
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Figure 4. Distribution dN/dD of distances in the 2MRS sam- 
ple compared to expectations for a homogeneous universe. The 
model curve is proportional to D 2 <f> sc \(D), where so i denotes 
the selection function. 



3 ACCURACY OF FLUX PREDICTIONS 

In this section we address fluctuations associated with 
the finite number of galaxies in flux maps modeled from 
a galaxy catalog. Although the sliding-box technique ef- 
ficiently suppresses these fluctuations, they still pose a 
potential limitation to the applicability of the method. 
This makes it important to quantify them, as we do in 
the following. 

First we briefly discuss the construction of flux maps 
from a galaxy catalog. We express the model flux from 
galaxy i as follows: 



Fi 



F WjJ{Di) 
4ttD 2 



(8) 



where Fo is a normalization constant, Wi is the weight as- 
signed to galaxy i by the sliding-box technique, Di is the 
galaxy distance, and J(D) represents the fraction of the 
integral flux from a source at distance D that survives 
attenuation by redshift and interaction with the ambi- 
ent matter. The function J(D) is different for UHECRs, 
neutrinos, and gamma rays; J(D) = 1 corresponds to 
no attenuation. To keep the discussion general we do not 
specify J(D) at this point. 

Generally speaking, the model flux for a given direc- 
tion on the sky is constructed by adding and averaging 
the fluxes of individual sources close to the line of sight. 
This can be done in various ways, in particular by dividing 
the sky into bins or by employing a smearing routine that 
distributes single-source fluxes over (part of) the sky. For 
the discussion of fluctuations the precise method is not 
very important; the critical parameter is the solid angle 
Af] over which flux contributions of individual sources 
are averaged. 

How large may fluctuations be in order not to spoil 
the accuracy of the flux map? The answer to this question 
depends clearly on the purpose of the map and can thus 
not be answered in general. A reasonable requirement, 
which we will pursue in the following, is that no significant 
contribution to the flux within the solid angle Afl should 
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come from a single source, i.e.0 



(9) 



where AF denotes the total flux of sources within Afl. 

With a few simplifying assumptions, the condition 
([9]) can be expressed in general terms. We first approx- 
imate AF ~ FAfi/n, where F = J^Fj stands for the 
total flux and f2 is the total solid angle occupied by the 
sample (fi = 4tt in the case of complete sky coverage). 
Equation ((9} then reduces to: 



An f 



(10) 



To satisfy this requirement it is necessary to have many 
sources in the solid angle AQ. However, this may be not 
sufficient because not all sources contribute the same flux 
at Earth. The situation is thus complicated by the depen- 
dence on distance D. Introducing the fraction f(D) of the 
total flux produced by the sources closer than D, and the 
number of these sources in the sample n(D), equation 
(|10[) can be rewritten in the following way, defining the 
quantity T: 



T 



AQ dn 

!Td/ >>L 



(ii) 



The number of sources n(D) is readily calculated from 
equation ©: 

r D 

n(D) = / dDD 2 v(D)L(M*(D)). 
Jo 

The fractional flux /(£>) = F(D)/F(F max ), where 

f r D 

F(D) = ^ J dDp{D)L(M,(D))w(D)J(D), 

represents the total flux from sources closer than D. 
Note the appearance of the weights w(D) assigned by the 
sliding-box method in the last equation. Neglecting devia- 
tions of v(D) from 1, which is reasonable on cosmological 
scales, we use the above equations to find that 



dn _ 3Ny(D) 
df ~ DJ(D) 



dDJ(D) (sliding box), (12) 



where Nv(D) = N(D, M*(D)) is the number of galaxies 
in the volume-limited sample at distance D (cf. eqs. lHJ). 
We stress that equation (I12|l is valid for fluxes modeled us- 
ing the sliding-box technique. For a single volume-limited 
subsample that is valid up to D max , a similar computation 
yields: 



dn 

w 



SN v (D n 



)D 2 



D?^J(D) 



dDJ(D) (vol. ltd.), (13) 



Comparing eqs (|12[) and (|13[1 . we see that the number of 
sources contributing to a given flux fraction at distance 
D is increased by a factor N v {D)D'i n ^/N v {D m ^)D A . 
This factor is unity at D — D m ax (where the sliding-box 
method offers no improvement), but may become very 
large at small distances. 

3 Alternatively, one could demand that Fj be smaller than the 
average flux per solid angle AC This criterion would lead to 
equation dlOI I directly. 



Inserting equation (|12|l into (|11[) brings us to the final 
expression for our criterion of small fluctuations: 



T 



3A» AV(F)rw Jo° max dPJ{D) 

n d d 

max 

J(D) 



> 1. 



(14) 



The three factors in equation (|14p respectively encode the 
dependence of T on the angular scale, on the statistics of 
the flux-limited sample, and on the attenuation of the 
model flux. The last factor reduces to unity in the case of 
no attenuation. The equation has to hold for all values of 
D; when it is violated 0(1) fluctuations in flux may occur 
in regions of angular size AO due to the contribution of 
a single source. 

We now discuss some of the quantities entering equa- 
tion (|14p . The number of galaxies Nv(D) becomes small 
at both very large and very small distances, potentially 
leading to large fluctuations. For large distances this can 
be prevented by considering only sources up to a maxi- 
mum distance -D ma x and assuming an isotropic flux from 
sources beyond that distance. Alternatively, particle hori- 
zons may provide a natural maximum distance (see be- 
low). In the case of small distances, the number of nearby 
sources is small while their contribution may be impor- 
tant due to their proximity. Unlike the fluctuations at 
large distances which are due to our poor knowledge of 
the galaxy distribution at those scales, the fluctuations 
at small distances are physical and may represent the ac- 
tual flux variations due to close sources. Their complete 
treatment may require a case-by-case study of the most 
nearby objects. 

The last factor in equation (|14p encodes the effect 
of flux attenuation, which can play an important role in 
modeling the flux of UHECRs and of very-high energy 
gamma rays. Focusing on the case of UHECRs, we show 
in figure [5] the flux attenuation factor J as a function of 
distance for UHECR protons. The attenuation factor is 
obtained us ing a numerical cosmic-ray pro pagation code 
described in lKoers fc Tinvakovl (|2008l . l2009h . For compar- 
ison the attenuation factor due to redshift only is also 
shown in the figure. In figure [6] we show the quantity 



J C °" dDJ(D) 
D msx J(D) 



(15) 



which accounts for flux attenuation in equation (|14[) . Note 
that, as indicated in figure[5] the horizon for UHECR pro- 
tons above 60 EeV is around 200 Mpc. Since sources be- 
yond this distance do not contribute to the observed flux, 
the requirement T>1 should be satisfied automatically. 
We observe from figure [6] that A indeed blows up around 
200 Mpc, which guarantees that T ^> 1 for any value of 
Afl or N v ■ 



4 EXAMPLE: UHECR FLUX PREDICTIONS 
USING THE 2MRS CATALOG 

In this section we apply the sliding-box technique to 
model the flux of UHECR protons with energies in excess 
of 60 EeV from sources tracing the distribution of matter 
in the Universe. In modeling the effect of flux suppression 
due to attenuation, we assume a power-law injection spec- 
trum with index p = 2.2 extending to very high energies. 
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Figure 5. Flux suppression factor J as a function of distance 
for three different scenarios: redshift only, UHECR protons 
with energy above 40 EeV, and UHECR protons with energy 
above 60 EeV. In producing this figure we have assumed a 
power-law injection spectrum with index p = 2.2 extending to 
very high energies. 
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Figure 6. Flux attenuation factor A defined in eq. H15H as a 
function of distance for the same scenarios as shown in figure [5] 



For comparison we also model the flux distribution using 
a single volume-limited galaxy sample up to 250 Mpc. 

The distribution of matter is modeled using a subset 
of the 2MRS galaxy sample. This subset does not cover 
the galactic plane with \b\ < 10°, nor the region with 
S < —30° (see footnote [T]). Due to these cuts the catalog 
covers 63% of the sky, so that the total field of view is O = 
0.63 ■ 4-7T = 7.9 srad. In figure [7] we show the number of 
galaxies in a volume- limited sample Ny(D) of the 2MRS 
as a function of distance D. 

In actual flux computations, dividing the sky into 
bins of fixed size has disadvantages related to boundary 
effects and the arbitrariness of the binning scheme. These 
problems are avoided with an angular smearing routine, 
which essentially replaces the point-source flux of an in- 
dividual source by a (Gaussian) probability distribution. 
Adopting a smearing routine, the flux in a given direction 
n is computed as follows: 



(16) 





50 



11)0 



200 



251) 



301) 



150 
D (Mpc) 

Figure 7. Size of a volume-limited sample of our 2MRS sample 
up to distance D. 



where 



F t e X p(-6 2 /8* 

7T<93 



(17) 



Here Fi represents the flux from galaxy i, 9 denotes the 
angle between the galaxy and the line of sight n, and f5 s 
is the smearing angle. 

In figure [8] we show model UHECR flux maps ob- 
tained with equation (|16p using the full 2MRS sample 
with the sliding-box method (top panel) and, for com- 
parison, using a volume-limited subsample extending to 
250 Mpc (bottom panel). A comparison between the two 
panels demonstrates the significant increase in accuracy 
achieved with the sliding-box technique. In particular, the 
shot noise artefacts that are visible in the bottom panel 
are absent in the top panel. 

We now consider the fluctuations associated with the 
finite number of galaxies in the 2MRS catalog. In figure [5] 
we show the quantity T defined in eq. (|14[) as a function of 
D for the exemplary case of 60 EeV UHECR protons. For 
comparison we also show T for the volume-limited sample 
(using eqs. (|11[) and (|13|l ). Because 63% of the total flux 
is contained within the opening angle 6 B , the solid angle 
AQ, that enters in eq. (|14p is related to 8 S as follows: 



AO, = 



2tt(1 



0.63 



(18) 



We observe from figure |5]that, for the sliding-box method, 
T > 1 for distance between 5 and 250 Mpc. Hence the 
results of the previous section imply that fluctuations as- 
sociated with the finiteness of the galaxy sample should 
be small. For the volume limited subsample, on the other 
hand, the figure indicates that strong fluctuations are to 
be expected. 

The estimates on the strength of fluctuations can be 
verified through a direct computation of the ratio of indi- 
vidual flux contributions (j>i to the total flux $. Sampling 
over many directions on the sky, we compute the total flux 
$ via equation (|16[) and keep track of c/wx = max (pi , the 
largest individual contribution to $. The distribution of 
</>max/ ( l > is shown in figure [TOl for the flux maps shown 
in figure |8j i.e. for the case of UHECR protons with en- 
ergy in excess of 60 EeV and smearing angle S = 3°. 
As may be verified from figure 1101 4>max/$ <C 1 for the 
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Figure 8. Aitoff projection of the sky in galactic coordinates 
showing the model flux of UHECR protons above 60 EeV from 
sources tracing the distribution of galaxies up to 250 Mpc. The 
grayscale shows the relative flux on a logarithmic scale. (Areas 
in white are not covered by our subsample of the 2MRS cata- 
log.) The top panel shows the flux constructed from the original 
flux-limited sample using the sliding-box technique; the bot- 
tom panel shows the flux constructed from a volume-limited 
subsample up to 250 Mpc. For both cases we have removed 
sources closer than 5 Mpc and smeared the flux distribution 
with 9 S = 3°. 
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Figure 9. T as a function of D for the case of 60 EeV UHECR 
protons modeled using our subset of the 2MRS catalog. The 
figure applies to a smearing angle 9 B =3°, corresponding to 
AQ = 0.01 srad; the total field of view f2 = 7.9 srad. 



sliding-box method: no single source outshines the bulk. 
On the other hand, if we model the flux distribution from 
the volume-limited sample (bottom panel of figure |5), we 
find that the distribution of <^ ma x/'& peaks near 1. In this 
case 0(1) fluctuations in the predicted flux due to a single 
source are common (which is also clear from the bottom 
panel of figure |SJ , which means that the galaxy sample is 
too small to provide an accurate flux map. 
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Figure 10. Distribution of F ma x/-FAf2, i.e. the ratio of the 
largest individual contribution to the total flux within solid 
angle AC This figure applies to the flux map shown in the top 
panel of figure [8] 



5 SUMMARY 

We have addressed the problem of flux calculation from 
sources tracing the galaxy distribution in the Universe. 
We have discussed a sliding-box weigh ting scheme, build- 
ing on the work o f lLvnden-Beil l|l97ll ). that makes use of 
the information contained in a flux-limited galaxy catalog 
in the most efficient way. This scheme allows us to rep- 
resent the distribution of matter up to large scales while 
representing the distribution of matter on small scales 
with maximum accuracy. The resulting weight function is 
related to the galaxy luminosity function and may be used 
to infer the latter under the assumption of its coordinate 
independence. 

The sliding-box weighting scheme suppresses effi- 
ciently fluctuations due to the finite size of the sample at 
most distances except the largest and the smallest ones. 
We have presented estimates on the size of the remaining 
fluctuations. These estimates can be used to determine a 
maximum distance at which the catalog should be cut, or 
to find the minimum angular scale on which flux maps 
can be constructed accurately. We would like to stress 
that our estimates regard the size of fluctuations, and not 
their importance. For example, in a statistical test based 
on model flux distributions, fluctuations of order unity 
may be acceptable if the overall flux distribution shows 
very strong contrasts or when the angular scales of inter- 
est are much larger than the scale at which fluctuations 
occur. 

An advantage of the sliding-box scheme is that it al- 
lows a straightforward generalization to the cases when 
the sources trace preferentially certain types of galax- 
ies, or when the source luminosity is correlated with the 
galaxy type. Such effects may be accounted for by pre- 
weighting the galaxies in the catalog in a corresponding 
way and modifying accordingly the sliding-box weighting 
scheme. 
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